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Abstract 

Current experimental data indicate that the Koide relation for the pole masses 
of charged leptons, which can be parametrized as C}f° le = 2/3, is valid up to 
the accuracy of 0(1O~ 5 ). We show that the running masses of charged leptons 
fail in satisfying the Koide relation (i.e., Qi{fj) ^ 2/3), but the discrepancy 
between Qi(fi) and Qf olc is only about 0.2% at fi = M z . The Koide-like 
relations for the running masses of neutrinos (1/3 < Q U {M Z ) < 0.6), up-type 
quarks (Qu(M z ) ~ 0.89) and down-type quarks (Q D (M Z ) ~ 0.74) are also 
examined from M z up to the typical seesaw scale M R ~ 10 14 GeV, and they 
are found to be nearly stable against radiative corrections. The approximate 
stability of QuO-0 anc ^ QdO-O * s niainly attributed to the strong mass hierarchy 
of quarks, while that of Q ; (/i) and Q u (fi) is essentially for the reason that the 
lepton mass ratios are rather insensitive to radiative corrections. 
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|~T] Thanks to the precise measurements of three charged lepton masses [1], 

m e = (0.510998918 ± 0.000000044) MeV , 
m M = (105.6583692 ± 0.0000094) MeV , 

m T = (l776.99l°;*j) MeV , (1) 

the empirical Koide mass relation [2] 

n poic _ m e + m„ + m T 2 

— -; ri — o y z > 

has been testified up to the accuracy of (9(10 -5 ). Namely, the above experimental data 
yield -0.00001 < Qf ole - 2/3 < +0.00002. This precision is so amazing that there might 
exist some underlying physics behind the Koide relation. A number of authors have tried 
to understand Eq. (2) and extend it to neutrino masses or quark masses based on possible 
flavor symmetries and phenomenological conjectures [3-6]. It is certainly impossible to get a 
universal Koide relation for both charged fermions and neutrinos by using the experimental 
data at low energy scales, nor is it likely to achieve such a mass relation at any high energy 
scales. In this category of exercises, one has to take care of the concepts of fermion masses 
and clarify the running mass behaviors from one energy scale to another. 

Note that Eq. (2) holds for the pole masses of three charged leptons, which have been 
denoted as m e , and m T . Is it also applicable to the running masses m e (/x), rn^ji) and 
m r (/i) at a given energy scale (e.g., \i = M z )l We are going to answer this question both 
analytically and numerically. 

If the Koide relation is not a universal relation for fermion masses at low energy scales, 
one may wonder whether a universal Koide-like relation is achievable at a superhigh energy 
scale. We shall show that the answer to this question is negative for both charged fermions 
and neutrinos, but the dynamical reasons are somehow different in these two cases. 

Let us define the following Koide-like parameters for the running masses of charged 
leptons (f — I), up-type quarks (/ = U), down- type quarks (/ = D) and neutrinos (f — v): 

QfW = r , == ,2 , (3) 



\/mJjj) + y/m y (fi) + y/m z (fj,) 



where 



(e, 


(i, 


r) , 


for / = 


-l; 


(1, 


2, 


3) , 


for / = 


v ; 


(u, 


c, 


t) , 


for / = 


U; 


(d, 


s, 


b) , 


for / = 


D . 



As pointed out by Gerard et al [6], Qj- = 1 (maximum) would hold when three fermion masses 
had the extremely hierarchical spectrum (e.g., m x m y m z ), while Qj- = 1/3 (minimum) 
would occur if three fermion masses were exactly degenerate (i.e., m x = m y = m z ). It might 
be a puzzle that Qf ole happens to take the mean value of QJ m and Q™ ax . The main purpose 
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of this paper is just to examine the deviations of realistic Q h Q u , Q v and Q D from Qf ole at 
and above the electroweak scale. 



First of all, let us calculate the discrepancy between Qt(M z ) and Qf olc . The running 
masses of three charged leptons are related with their corresponding pole masses m l 

in a very simple way [8] , 

m M = m A l -^i) , (5) 

where 



a(n) 



7T 



- in — — — + 1 

2 m, (^) 



(6) 



with a(/i) being the fine-structure constant at the energy scale ji. It is then straightforward 
to establish the relationship between Qi(fi) and Qf ol ° by combining Eqs. (3) and (5). Taking 
account of the strong mass hierarchy m e <C <C m T (more explicitly, m e /m ~ 0.00484 
and m^/m T m 0.0595 [1]), we find that it is more instructive to make the following analytical 
approximation: 



(7) 



One can clearly see that the deviation of Qi(/u) from Qf olc is strongly suppressed, due to the 



smallness of ^Jm^/m T and that of A,. 

At n = M z with a(M z ) = (128.89) -1 [9], the running masses of three charged leptons 
can be evaluated by solving Eq. (5) with the inputs of the pole masses given in Eq. (1). 
For our purpose, we only need to make use of the central values of m e , and m T in the 
calculation. The results are 

m e (M z ) = 0.486755106 MeV , 
m^Mz) = 102.740394 MeV , 

m T (M z ) = 1746.56 MeV . (8) 

With the help of Eq. (6), we obtain = 0.0276 and A r = 0.0171 at M z . Therefore, 

Q t (M z ) = 1.00188 x Qf olc = 0.66792 , (9) 

achieved from Eqs. (3) and (5). If the analytical approximation made in Eq. (7) is used to 
estimate Q t (M z ), one can get Q t (M z ) ~ 1.00256 xQf olc = 0.66837. This result is compatible 
with Eq. (9) and implies that Eq. (7) is actually a reasonable approximation. 

Obviously, the discrepancy between Q t (M z ) and Qf olc is only about 0.2%. This tiny 
difference makes sense in physics, because it is much larger than the accuracy of Eq. (2). 
In other words, the possibility of Q t (M z ) = 2/3, which inversely leads to Qf olc — 2/3 = 
—0.00188, has been ruled out by the experimental data. We are therefore left with the 
conclusion that the running and pole masses of three charged leptons cannot simultaneously 
satisfy the Koide relation. 
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Now let us turn to the Koide-like relations of quark masses. We concentrate on the 



running masses of six quarks at /i — M z [10] 1 , 



m u {M z ) = (1.7 ±0.4) MeV , 
m c (M z ) = (0.62 ± 0.03) GeV , 

m t {M z ) = (171 ± 3) GeV ; (10) 



and 



m d {M z ) = (3.0 ± 0.6) MeV , 
m s (M z ) = (54 ± 11) MeV , 

m b (M z ) = (2.87 ±0.03) GeV . (11) 

The Koide-like parameters Q v and Q D defined in Eq. (3) can then be calculated at M z 
with the help of Eqs. (10) and (11). For simplicity, only the central values of those quark 
masses are taken into account. The results are 

Q V {M Z ) « 0.89 , Q D (M Z ) « 0.74 . (12) 

We see that both numbers significantly deviate from 2/3 at the electroweak scale. 

To examine whether the Koide-like relation Q V (M X ) as Q D (M X ) 2/3 could hold at a 
superhigh energy scale M x (e.g., M x ~ X0 14-16 GeV), one may make use of the one-loop 
renormalization-group equations (RGEs) of quark Yukawa couplings [12] in the standard 
model (SM) or in the minimal supersymmetric standard model (MSSM). Note that the 
strong hierarchy of charged fermion masses and that of the quark mixing angles allow us 
to simplify those RGEs to a great extent [13]. In particular, the RGE evolution of m u /m c , 
m d /m s and m e /m ll from M z to M x are negligibly small in both the SM and the MSSM. 
Radiative corrections to the mass ratios m c /m t , m s /m b and m fl /m T can be written as 



™ C (M X ) _ m c (M_ 



z) 



m t (M x ) m t (M z ) 
m s (M x ) _ m s (M z ] 



m b (M x ) m b (M z ) 
m^(M x ) m^M z ] 



m T (M x ) m T (M ZJ 
where the evolution functions Xu> Xd anc ^ Xi are defined by 



Xi , (13) 



X/ = exp 



(14) 



We do not consider the pole masses of six quarks, because they are neither directly measurable 
nor relevant for model building. In particular, the pole masses of three light quarks (u, d, s) involve 
large uncertainties, since they can only be evaluated in the region with a large a s (fi) [11]. 
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with t = ln(/x), / = (U,D,Z) and y i 



t, b, t) being the Yukawa couplings. In the SM, 



b D = c t = —3/2 and c v = c r 



have a T 



-3, W 



-1 and c 



u 



6, = hold; while in the MSSM, we 
a l = b l = 0. The running behaviors 



of Xf are typically illustrated in Fig. I. One can see that Xi is insensitive to radiative 
corrections in both the SM and the MSSM. In contrast, Xu ma Y significantly depart from 1 
when M x M z holds, especially for sizable tan/3 in the MSSM. Note that Xd > 1 holds 
in the SM, as a consequence of a D > in Eq. (14). 

With the help of Xf (f° r / — U,D,Z) given above, it is easy to derive the Koide-like 
parameter Qf{M x ) in terms of Qj(M z ). We approximately obtain 



Qv(M x ) 


~ Qv( M z) 


1 + 2 (l - 




Qb(M x ) 


~ Qb(M z ) 


l + 2(l- 




Qi(M x ) 


~ g z (M z ) 


l + 2(l- A 


/*)\ 



m c (M z ) 
m t (M z ) 



m s (M z 
\ m b {M z 



m T (M z ) 



(15) 



One can see that the deviation of Qj(M x ) from Q f (M z ) is significantly suppressed, because 
of (i) the strong mass hierarchies of charged fermions and (ii) the small departure of Xf 

from 1. Indeed, ^m c /m t 0.06, ^Jm s /m b ps 0.14 and y / m /J /m T ss 0.24 at M z . Taking 
M x = M R ~ 10 14 GeV (the typical scale of heavy right-handed Majorana neutrinos in the 
seesaw models [14]), for example, we get y/x~v ~ 0.97, t/Xd ~ 1-03 and ^fxx ~ 1-00 in the 

w 0.98 in the MSSM with tan/3 = 50. Therefore, 



SM; or ^ ~ 0.86, ^ ~ 0.91 and 
Q f (M R )-Q f (M z ) 

in the SM; or 2 

Q f (Mj-Q f (M z ) 



3.3 x 10" 3 , 
-6.5 x 10~ 3 
0, 



1.4 x lO- 2 
1.8 x 10~ 2 
4.8 x 10~ 3 



for / = U , 
for / = D , 
for / = I 



for / = U , 
for / = D , 
for / = I 



(16) 



(17) 



in the MSSM with tan/3 = 50. It is obvious that radiative corrections to the Koide-like 
parameters are very small. 



Finally, we consider the Koide-like relation in the neutrino sector. Although the pole 
masses of three neutrinos are different from their running masses at fj, — M z , this difference 



2 The deviation of Qj(M R ) from Qj(M z ) in the MSSM will become much milder, if tan/3 takes 
smaller values (tan/3 = 10 shown in Fig. 1(b), for example). 
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is negligibly tiny because it is strongly suppressed by the Fermi coupling constant. For 
simplicity, we mainly calculate the Koide-like parameter Q U (M Z ) and examine its sensitivity 
to radiative corrections from M z up to the seesaw scale M R . 

Note that the absolute values of three neutrino masses remain unknown, but their upper 
bound is expected to be m i < 0.23 eV (for % = 1,2,3) [15]. A global analysis of current 
experimental data on neutrino oscillations yields [16]: Am^ = m 2 — m\ = (8.0 ± 0.3) x 
10~ 5 eV 2 and Am| 2 =m\-m\ = ±(2.5 ±0.3) x 10~ 3 eV 2 at the 99% confidence level. Since 
the sign of Am§ 2 has not been fixed, the neutrino mass spectrum can be classified into two 
general categories: 

• Normal hierarchy: m 1 < m 2 < m 3 , including the possibility that three neutrino masses 
are nearly degenerate (Am| 2 > 0). 

• Inverted hierarchy: m 3 < m 1 < m 2 , including the possibility that three neutrino 
masses are nearly degenerate (Am 32 < 0). 

Allowing m 1 or m 3 to vary up to its upper bound and inputting the experimental values 
of Am 21 and Am| 2 , we may use Eq. (3) to calculate Q U (M Z ). Our numerical results are 
illustrated in Fig. 2. One can see that the upper limit of Q U (M Z ) is achieved at m 1 = for 
the normal neutrino mass hierarchy or at m 3 = for the inverted neutrino mass hierarchy. 
Indeed, 



Qu(M z ) 



^Am| n 


/Am 2 21 + 


Amis | 


^Am^ + 1 


/Am 2 ! + 


|Am§ 2 | 


^\Am 2 32 \ 


\/l Am 32| 


- Amli 



(18) 



1 ± 2 V = , (m 3 -> 0) . 



Am 32 | ± w|Am 32 — Am 21 



It is then straightforward to understand Q U (M Z ) ~ 0.6 for m 1 — > and Q U (M Z ) m 0.5 for 
m 3 — > (as shown in Fig. 2) in the approximation of |Am 32 | ^> Am 21 . When three neutrino 
masses are nearly degenerate, Q V (M Z ) approaches its minimal value 1/3 no matter whether 
the sign of Am| 2 is positive or negative. 

The running masses of three neutrinos at M R can be evaluated by using the one-loop 
RGEs given in Ref. [12]. In order to understand the sensitivity of Q U (M R ) to radiative 
corrections, it is more convenient to consider the RGE running behaviors of m i jm^ (for 
i j)- We obtain 

m {M x ) ~ m 3 {M z )^ 1 ) 

in the approximation of r-lepton dominance [17], where 3 



3 This formula is valid for Majorana neutrinos. If neutrinos were Dirac particles, the coefficient 
C/(8vr 2 ) on the right-hand side of Eq. (20) should be replaced by C/(16vr 2 ) [17]. 
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X ( „ 3) = exp 



(J AnM 



8tt 2 



f M y 2 r l^l 2 -l^l 2 )dt 



(20) 



with C = -1.5 in the SM or C = 1 in the MSSM and U 3i (for i = 1,2,3) being the 
elements of the 3x3 lepton flavor mixing matrix. Because of y^/(8ir 2 ) ~ 1.3 x 1CT 6 (SM) or 
^/(8tt 2 ) « 1.3x 10~ 6 (1 + tan 2 /3) (MSSM) at M z , together with \U 3i \ 2 < 1, the departure of 
from 1 is expected to be tiny. This observation is illustrated in Fig. 3, where 9 U ~ 33.8° 
and 9 ~ 45° [17] have typically been input for the matrix elements \U 31 \ = sin 9 V sin 9, 



\U 32 \ = cos 9 U sin 9 and |C/ 33 | = cos 6*. 



Note that ~ 1 is essentially independent of the possible mass hierarchies of three 
neutrinos [18]. It turns out that Q U (M X ) fa Q U (M Z ) is a very good approximation. In view 
of Fig. 2, we conclude that there is no hope to achieve the Koide-like relation Q u (fi) ~ 2/3 
for neutrino masses. 



In summary, the updated Koide relation for the pole masses of charged leptons 
satisfies Qf ole = 2/3 at the precision level of O(10~ 5 ). This amazing accuracy motivates 
us to examine whether the running masses of charged fermions and neutrinos have the 
similar Koide relation at a given energy scale. We have shown that the running masses of 
charged leptons cannot satisfy the Koide relation (i.e., Qi(fi) ^ 2/3), but its discrepancy 
from Qf ole is only about 0.2% at fi = M z . The Koide-like relations for the running masses 
of neutrinos (1/3 < Q U (M Z ) < 0.6), up-type quarks (Q V (M Z ) ~ 0.89) and down-type 
quarks (Q D {M Z ) ~ 0.74) are analyzed from M z up to the typical seesaw scale M R ~ 10 14 
GeV. We find that they are nearly stable against radiative corrections, just like Q ; (/x). 
The approximate stability of Q v (fi) and Q D (/i) are mainly attributed to the strong mass 
hierarchy of quarks. In contrast, the approximate stability of Q;(/i) and Q„(/i) is essentially 
for the reason that the lepton mass ratios are very insensitive to radiative corrections. 

Although it is impossible to get a universal Koide relation for both charged fermions and 
neutrinos at a given energy scale, our work remains useful for model building in order to 
explore the underlying similarities and differences between lepton and quark masses. For 
example, 

Q V (M Z ) > Q D (M Z ) > Q,(M Z ) > Q U (M Z ) (21) 

is an interesting result of our numerical analysis. Once the absolute scale of neutrino masses 
is measured and the value of Q V (M Z ) is fixed, one may then speculate whether these four 
parameters could be related with one another in a phenomenological way. On the other hand, 
we remark that more theoretical efforts are needed to look for possible flavor symmetries 
and their breaking effects behind the Koide-like relations. 
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FIG. 1. Illustration of Xu> Xd an d Xi changing with the energy scale from M z to M R . 
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FIG. 2. The Koide-like parameter Q U (M Z ) evaluated by using current neutrino oscillation data. 
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FIG. 3. Illustration of xh changing with the energy scale from M z to M R . The solid and 
dashed lines denote the Majorana and Dirac cases, respectively. We have typically input tan (3 = 50 
and m l = 0.2 eV with m\ 2 > in our calculation. The almost identical result can be obtained 
when m\ 2 < is taken. Note that the deviation of from 1 will be much milder, if tan f3 takes 
smaller values. 
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